Surface magnetoinductive breathers in two-dimensional magnetic metamaterials 
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We study discrete surface breathers in two-dimensional lattices of inductively-coupled split-ring 
resonators with capacitive nonlinearity. We consider both Hamiltonian and dissipative systems and 
analyze the properties of the modes localized in space and periodic in time (discrete breathers) lo- 
cated in the corners and at the edges of the lattice. We find that surface breathers in the Hamiltonian 
systems have lower energy than their bulk counterparts, and they are generally more stable. 
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PACS numbers: 63.20.Pw, 75.30.Kz, 78.20.Ci 

Theoretical results on the existence of novel types of 
discrete surface solitons localized in the corners or at 
the edges of two-dimensional photonic lattices [J H, [1[ 
have been recently confirmed by the experimental obser- 
vation of two-dimensional surface solitons in optically- 
induced photonic lattices [J] and two-dimensional waveg- 
uide arrays laser- written in fused silica [a, [y] . These two- 
dimensional nonlinear surface modes demonstrate novel 
features in comparison with their counterparts in trun- 
cated one-dimensional waveguide arrays |7|, |8|, |9| . In par- 
ticular, in a sharp contrast to one-dimensional surface 
solitons, the mode threshold is lower at the surface than 
in a bulk making the mode excitation easier [2[ . 

Recently, it was shown |l_0(j that, similar to discrete 
solitons analyzed extensively for optical systems, surface 
discrete breathers can be excited near the edge of a one- 
dimensional metamaterial created by a truncated array 
of nonlinear split-ring resonators. Networks of split-ring 
resonators (SRRs) that have nonlinear capacitive ele- 
ments can support nonlinear localized modes or discrete 
breathers (DB's) under rather general conditions that de- 
pend primarily on the inductive coupling between SRRs 
and their resonant frequency (TH fl21 | . The correspond- 
ing one-dimensional surface modes have somewhat lower 
energy (in the Hamiltonian case) and can easily be gen- 
erated in one-dimensional SRR lattices [lfj ■ 

In this Brief Communication, we develop further those 
ideas and analyze two-dimensional lattices of split-ring 
resonators. Similar to the optical systems, we find that 
two-dimensional lattices of inductively-coupled split-ring 
resonators with capacitive nonlinearity can support the 
existence of long-lived two-dimensional discrete breathers 
localized in the corners or at the edge of the lattice. 

We consider a two-dimensional lattice of SRRs in both 
planar and planar-axial configuration [see Figs. QJa,b)]. 
In the planar configuration, all SRR loops are in the same 



plane with their centers forming an orthogonal lattice, 
while in the planar-axial configuration the loops have a 
planar arrangement in one direction and an axial config- 
uration in the other direction. Each SRR is equivalent 
to a nonlinear RLC circuit, with an ohmic resistance R, 
self-inductance L, and capacitance C. We assume that 
the capacitor C contains a nonlinear Kerr-type dielectric, 
so that the permittivity e can be presented in the form, 



e(|E| 2 ) =e [e e + a 
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where E is the electric field with the characteristic value 
E c , ei is linear permittivity, eo is the permittivity of 
the vacuum, and a = +1 (—1) corresponding to self- 
focusing (self-defocusing) nonlinearity, respectively. As 
a result, each SRR acquires the field-dependent capac- 
itance C(|E| 2 ) = e(|E g | 2 ) A/dg, where A is the area of 



the cross-section of the SRR wire, 



induced along the SRR slit, and d g 



E ff is the electric field 



is the size of the slit. 



The field E g is induced by the magnetic and/or electric 
component of the applied electromagnetic field, depend- 
ing on the relative orientation of the field with respect to 
the SRR plane and the slits [13| . Below we assume that 
the magnetic component of the incident (applied) electro- 
magnetic field is always perpendicular to the SRR plane, 
so that the electric field component is transverse to the 
slit. With this assumption, only the magnetic compo- 
nent of the field excites an electromotive force in SRRs, 
resulting in an oscillating current in each SRR loop. This 
results in the development of an oscillating voltage differ- 
ence U across the slits or, equivalently, of an oscillating 
electric field E g in the slits. 

If Q is a charge stored in teach SRR capacitor, from 
a general relation of a voltage-dependent capacitance 
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(a) 



(b) 



FIG. 1: Schematic of a two-dimensional lattice of split-ring 
resonators for (a) planar and (b) planar-axial geometries. In 
both the geometries the magnetic component of the applied 
field is directed along the SRR axes, while the electric field 
component is transverse to the slits. 



C(U) = dQ/dU and Eq. ((TJ), we obtain 

U 2 



= C e 1 



(2) 



where U — d g E g , Cg = eoet(A/d g ) is the linear capaci- 
tance, and U c = d g E c . We assume that the arrays are 
placed in a time- varying and spatially uniform magnetic 
field of the form 



H = Ho cos(wt), 



(3) 



where Ho is the field amplitude, u) is the field frequency, 
and t is the time variable. The excited electromotive 
force £ , which is the same in all SRRs, is given by the 
expression 



£ = So sin(wi), £q = fio lo S Hq, 



(4) 



where S is the area of each SRR loop, and no the per- 
mittivity of the vacuum. Each SRR exposed to the ex- 
ternal field given by Eq. © is a nonlinear oscillator 
exhibiting a resonant magnetic response at a particular 
frequency which is very close to its linear resonance fre- 
quency cue = 1 / y/L Ce (for R ~ 0). 

All SRRs in an array are coupled together due to mag- 
netic dipole-dipole interaction through their mutual in- 
ductances. However, we assume below only the nearest- 
neighbor interactions, so that neighboring SRRs are cou- 
pled through their mutual inductances M x and M y . This 
is a good approximation in the planar configurations [see 
FigQJa)], even if SRRs are located very close. Validity of 
the nearest-neighbor approximation for the planar-axial 
configuration [see FigQTb)] has been verified by taking 
into account the interaction of SRRs with their four near- 
est neighbors. Assumimg that the mutual inductance 

(s) 

Mx.y between an SRR and its s— th neighbor decays with 
distance as M^l — M x ^ y /s 3 [l2j|, we find practically the 
same results. Therefore, the electrical equivalent of an 
SRR array in an alternating magnetic field is an array 
of nonlinear RLC oscillators coupled with their nearest 
neighbors through their mutual inductances; the latter 
are being driven by identical alternating voltage sources. 



Equations describing the dynamics of the charge Q n ,m 
and the current I n . m circulating in the n, m— th SRR 
may be derived from Kirchhoff's voltage law for each 
SRR [UG1 
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where f(Q n ,m) = Un.m is given implicitly from Eq. f2}. 
Using the relations 

uj^ 2 = LCt, t = tuje, I c = U c Lu e Ce, Q c = Ce~U c (7) 
and Eq. Q, we normalize Eqs. © and ^ to the form 

dQn,m 

dr 

din,m ■ , ft \ i \ (din—lm di n _i_i r 

+ J In.m + j(.Qn,m) + K 



(9) 
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£ sin(Or), (10) 



where 7 = RCguj£ is the loss coefficient, X XjV = 
M XtV /L are the the coupling parameters in the x— and 
y— direction, respectively, and £0 = £o/U c . Note that the 
loss coefficient 7, which is usually small (7 <C 1), may 
account both for Ohmic and radiative losses [14]. Ne- 
glecting losses and without applied field, Eqs. © and 
(JTU]) can be derived from the Hamiltonian 



m T *n,m 



n.m ^ 

^ ^ {Aa; Qn,m Qn-\-l,rn Ay Qn,m Qn,m-\-l }, (11) 

where the nonlinear on-site potential V n ,m is given by 

V4,m = V{q n , m ) = j f(<l'n,m) d q'n, m , ( U ) 



and q n . rn = dq nm /dr . Analytical solution of Eq. © for 
Un,m = /(^n,m) with the conditions of Un,m 

being real 

and u n m 

(q n ,m = 0) = 0, gives the approximate expres- 
sion 



f(q n ,m) ^ qn,m ~ ^?n,m + 3 {j^) q ™,™> ^ 

which is valid for relatively low q n (q n < 1, n = 
1, 2, N). Thus, the on-site potential is soft for a = +1 



and hard for a = — 1. In the 2D case the mutual induc- 
tances M x and M y may differ both in their sign, depend- 
ing on the configuration, and their magnitude. Actually, 
even in the planar 2D configuration with d x = d y a small 
anisotropy should be expected because we consider SRRs 
having only one slit. This anisotropy can be effectively 
taken into account by considering slightly different cou- 
pling parameters X x and X y . The coupling parameters 
\ x , y as well as the loss coefficient 7 can be calculated 
numerically for this specific model with high accuracy. 
However, for our purposes, it is sufficient to estimate 
these parameters for realistic (experimental) array pa- 
rameters, ignoring the nonlincarity of the SRRs and the 
effects due to the weak coupling as in Refs. [U, [T3| with 
the following typical values A as 0.02 and 7 w 0.01 . 

We construct discrete breathers located in the corner 
of a two-dimensional lattice of 15 x 15 sites using the 
anti-continuous limit method as in Ref. [ll| for the set 
of Eqs. (|9^- (fl"0|) . setting 7 = and eq — (Hamiltonian 
discrete breathers). For the case of a = +1 correspond- 
ing to self-focusing nonlinearity and period Tb = 6.69, 
we may construct linearly stable breathers for param- 
eters up to \ x = X y = 0.029. Breather stability has 
been checked through the Floquet monodromy matrix 
throughout the paper. For the case where an anisotropy 
is introduced, A^ < A y , linearly stable discrete breathers 
can be constructed up to A^ = 0.028 and simultaneously 
A y = 0.031, or for the case of planar-axial configuration 
up to X x = 0.031 and at the same time X y = —0.028. If 
we look for discrete breathers constructed in the middle 
of the upper edge of the lattice for example, we find that 
the values of the coupling where an instability occurs are 
slightly decreased (e.g. the upper stability limit of cou- 
pling for planar geometry is X x = X y = 0.028). Several 
cases of linearly stable discrete breathers are shown in 
Fig. [2] for a = +1. The same analysis holds for a = — 1 
(defocusing nonlincarity) where the upper stability limit 
for the values of couplings are of the same order of mag- 
nitude as for a = +1, both for the corner and edge 
breathers (see Fig. [3]) . The breather period in the lat- 
ter case is Tb = 5.8. 

Localized modes in the damped-driven case are con- 
structed for 7 = 0.01, £0 = 0.04 and a = +1 with 
the method described in Ref. [ill- The resulting lo- 
calized modes are called dissipative breathers, and their 
examples are shown in Fig. 0] for 7], = 5.8 and (a) 
X x = X y = 0.0007, (b) A^ = 0.0022 and X y = 0.0052, 
and (c) X x = 0.0052 and X y = -0.0022. The dissipative 
modes have been evolved in time, and we found that at 
long times some dissipative breathers constructed for rel- 
atively large couplings loose their initial shape and finally 
decay. 

Additionally, we calculate the total energy of discrete 
breathers in a lattice with planar and planar-axial con- 
figuration for a = +1 and T& = 6.69 (Hamiltonian case). 
Figure [5] shows the energy histograms of the relevant 
corner of the lattice normalized to the energy of the 
corner (1,1) breather. In order to construct the his- 




FIG. 2: Density amplitudes q n ,m for discrete Hamiltonian 
breathers constructed in (a-c) upper left corner or (d-f) upper 
edge of the lattice of 15 x 15 sites, a = +1 and Tj, = 6.69. 
(a,c) X x = X y = 0.028, (b,e) A* = 0.026 and \ v = 0.029, 
(c,f) X x = 0.029 and X y = -0.026. All plots depict a 5 x 5 
sublattice that includes the breather zones. 
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FIG. 3: Density amplitude q n , m for discrete Hamiltonian 
breathers constructed in (a-c) upper left corner or in (d-f) 
upper edge of a lattice of 15 x 15 sites, a = — 1 and Tb = 5.8. 
(a,d) X x = X y — 0.030, (b,e) A* 



(c,f) X x 



0.028 and A„ 



-0.025. All plots depict the 5 x 5 



sublattice around the linearly stable breathers. 
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FIG. 4: Density amplitude q n , m for discrete dissipative 
breathers for 7 = 0.01, e = 0.04, a = +1 and T b = 6.82, 
constructed in the upper left corner for (a) X x = X y = 0.0007, 
(b) A^ = 0.0022 and X y = 0.0052, and (c) X x = 0.0052 and 
A y = —0.0022. All plots depict the 5x5 sublattice around 
the breather. Dissipative breathers are very narrow and es- 
sentially confined on one lattice site. 



tograms centered in each of the lattice sites, we normal- 
ized it to the edge breather energy. In the case (a) the 
discrete breather is constructed in a lattice of coupling 
A x = Ay = 0.028, in (b) the case with anisotropy in cou- 
plings X x = 0.026 and X y = 0.029, while in the case (c), 
couplings are A^ = 0.029 and X y = —0.026. The energy 
of the discrete breathers as a function of the lattice site 
increases, i.e, as the discrete breather is constructed in 
the interior of the lattice energy is larger compared to 
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the discrete breather that is located in the corner of the 
lattice. 




FIG. 5: Histogram of the breather Hamiltonian. Breather 
difference energies AE for a — +1, Tb — 6.69 constructed 
in the upper left 3x3 corner of the lattice. Case (a) A^ = 
X y = 0.028, case (b) X x = 0.026 and X y = 0.029, and case (c) 
A^ = 0.029 and X y — —0.026. To evaluate AE we calculate 
the energy of the breathers centered at different sites and 
subtract the energy of the corner breather. 




FIG. 6: Amplitudes q n , m of the breather for a — +1, Tb = 
6.69 and X x — X y — 0.028, constructed on the site (1,1) for 
(a) t=0 and (b) t = 1450Tb, and the breather constructed on 
site (3,3), for (c) t = and (d) t = 1450T 6 . 

We note that in the one-dimensional case the bulk 
breathers have lower energy compared to the surface 



ones [lOj while in two-dimensional lattice the behavior 
is the contrary. We thus find that two-dimensional sur- 
face and especially edge breathers form easier. 

Finally, we study the time evolution of the discrete 
breather that is constructed in the corner site (1,1) and 
compare this case with a discrete breather centered at the 
(3,3) site for the coupling = X y = 0.028. The breather 
of the latter case after t — 95Xf, starts to loose its shape, 
in contrast to the breather of (1,1) site which survives 
for much longer times, viz. t = 1450Tfc [see Fig. 
For different coupling values such as X x = X y = 0.01 we 
find that both the corner (1,1) and inner (3,3) breathers 
remain stable for at least t = 1450T;,. This feature, while 
compatible with the fact that the corner breathers are 
more stable than inner ones, shows additionally that in 
finite lattices small changes in parameters may affect the 
stability properties of the breathers [l5[ . 

In conclusion, we have studied surface discrete 
breathers located in the corner and at the edge of the 
two-dimensional lattices of the split-ring resonators. Us- 
ing standard numerical methods, we have found nonlin- 
ear localized modes both in the Hamiltonian and dissipa- 
tive systems. Two-dimensional breathers in conservative 
lattices have been found to be linearly stable for up to 
certain (large) values of the coupling coefficient, in both 
planar and planar-axial configurations of the split-ring- 
resonator lattices. Dissipative discrete surface breather 
can retain their shapes for several periods of time, and 
they depending critically on the lattice coupling. Finally, 
we have found that the discrete breathers located deep 
inside the lattice have higher energy compared to the 
breathers located in the corners and at the edges. This 
distinct two-dimensional feature of nonlinear localized 
modes contrasts with the one-dimensional behavior being 
attributed to the larger number of neighbors of the two- 
dimensional lattice. Furthermore, the two-dimensional 
breathers located inside the lattice loose rapidly their 
initial shape as they evolve in time while the surface 
breathers are seen to be stable at least for t w 1500T;,. 
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